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Abstract 

A theory of the electronic structure and excitonic absorption spectra of PbS and PbSe nanowires 
and nanorods in the framework of a four-band effective mass model is presented. Calculations con- 
ducted for PbSe show that dielectric contrast dramatically strengthens the exciton binding in 
narrow nanowires and nanorods. However, the self-interaction energies of the electron and hole 
nearly cancel the Coulomb binding, and as a result the optical absorption spectra are practically 
unaffected by the strong dielectric contrast between PbSe and the surrounding medium. Mea- 
surements of the size-dependent absorption spectra of colloidal PbSe nanorods are also presented. 
Using room-temperature energy-band parameters extracted from the optical spectra of spherical 
PbSe nanocrystals, the theory provides good quantitative agreement with the measured spectra. 
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I. INTRODUCTION 



Solution-based chemical synthesis of semiconductor nanostructures has allowed tremen- 
dous flexibility in crystal morphology. After much work on zero-dimensional (OD) nanocrys- 
tals (NCs), attention is shifting to one-dimensional (ID) nanorods (NRs) and nanowires 
(NWs)i~- and the variation of material properties in the transition from OD to ID. The 
electronic structure of these crystals is the foundation for understanding their properties. 
Previously, the electronic structure of ID nanocrystals has been modeled using a variety of 
methods, including effective-mass theories based on k ■ p Hamiltonians^-, pseudopotential 
techniques^— , tight binding models^r— , and density functional theory*^ - — . The relax- 
ation of confinement in going from OD to ID goes hand-in-hand with an increase in the 
importance of Coulomb effects mediated through the nanocrystal's dielectric environment^. 

Lead-salt (PbS, PbSe, PbTe) nanocrystals offer unique advantages to study the interplay 
of these two effects. Their large exciton Bohr radii places them at the limit of strong 
confinement, while their large dielectric constants coupled with their mirror-like electron and 
hole spectra substantially reduce the Coulomb interaction in spherical quantum dots^ 1 ^. 
However, in a ID structure the Coulomb interaction can act primarily through the host 
medium, so it will not be screened as effectively as in OD^. Thus, the lead salts provide a 
unique system to study the transition from strong confinement to strong Coulomb binding 
as the length of the nanocrystal changes. 



NWs and NRs sur- 
7|. A type of adi- 



Within kp theory, the general treatment of the optical properties o: 
rounded by media with small dielectric constant was developed in Refs. 
abatic approximation naturally separates the calculation into pieces. In recognition of strong 
confinement perpendicular to the NR or NW axis, one first calculates the ID subband ener- 
gies and wavefunctions, while neglecting the Coulomb interaction. Next, using these wave 
functions of transverse electron and hole motion, one can calculate the longitudinal motion 
of the exciton, including corrections from image forces in the surrounding medium. To do 
that, the three-dimensional Coulomb potential is averaged to a one-dimensional Coulomb 
interaction between the electron and hole along the NW or NR axis. Using this potential, 
the spectra of ID excitons and their transition oscillator strengths are found. Finally, in 
NRs one should find the spectrum of the exciton center of mass motion, in order to include 
this additional effect of confinement. The main aspects of this framework were performed 
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for lead-salt nanowires recently by Rupasov^, although approximations to the simplified 
band structure used in that paper preclude the description of real experimental results. 

In this paper we present calculations of the ID subband energy spectra of lead-salt 
nanowires with arbitrary axis orientation, taking into account the multi- valley structure and 
accurate electron and hole energy-level dispersions in these semiconductors. For PbSe NWs 
with axis along the (100) direction, we calculate the spectra of ID excitons including self- 
interaction corrections. Surprisingly, the calculations show that although the binding energy 
of excitons in the smallest NWs reaches 350 meV, the optical transition energies are not 
affected by the small dielectric constant of the surrounding medium and are almost identical 
to the transitions between non-interacting electron and hole subbands. The cancelation 
of the exciton binding energy and the self-interaction corrections to the electron and hole 
levels is a consequence of the almost mirror symmetry of the conduction and valence bands 
of PbSe. The theoretical results agree well with the measured absorption spectra of (100) 
PbSe NRs. 

The paper is organized as follows. In Section [TT] we will describe the Hamiltonian governing 
the ID nanowire system, with solutions in Section IIHI In Section |IV] we present the effects 
of dielectric confinement and Coulomb forces on the ID exciton, with ID wavefunction 
solutions in Section [V] Experimental data and comparison with theory are presented in 
Section Wl\ followed by discussion and conclusion. 



II. FOUR BAND EFFECTIVE MASS MODEL 



PbS, PbSe, and PbTe are direct-gap semiconductors, with extrema of the conduction and 
valence bands at the L points in the Brillouin zone. The energy bands near the L point 
can be well-described within the four-band k • p model^ 1 ^. This model takes into account 
the direct interaction between the nearest conduction and valence bands, as well as the 
contributions of the remote bands to the electron and hole effective masses. Following Ref. 



24| , we use the multiband effective mass approximation and expand the full wave functions 



inside the nanorod as 



*(r)= E E WW> (!) 

Ai=±l/2 fi=±l/2 
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where \Lq ) and |Lg „) are the Bloch functions of the conduction band and valence band 
edge, respectively, at the L-point. The upper sign " ± " in the notation reflects the invariance 
of these functions with respect to the operation of spatial inversion. The smooth functions 
ty±i/ 2 { r ) an d ^± 1 / 2 ( r ) are the components of the conduction band and valence band 
spinor envelopes, respectively: 

/ \T/ C \ / \V" \ 

(2) 





' *S/2 ' 















The bi-spinor envelope function \1/ = I I is the solution of the Schrodinger equation 
H[p)^ = E^> . where p = hk = —ih'V is the momentum operator, and the Hamiltonian 



H(P) of Ref. 24| can be written in compact form as 



H(p) 



( (Eg Pi Pi \ f Pi _ Pt . , \ 

2 2m, 2m t J m mo 



Pi ^ Pt r * , (Eg Pi P\ \ f 

\ m m \ 2 2m, 2m^ / / 



(3) 



In Eq. ([3]) C/2 is the 2x2 unit matrix, a = \o X i^y^z\ are the Pauli matrices that 
act on the spinor components of the wave functions (/i = ±1/2), E g is the bulk energy 
gap, E is the electron or hole energy measured from the middle of the gap, mo is the 
free electron mass, p\ = Pi + Py , (P±d'±) = p x a x + pycr y , P t and Pi are the transverse 
and longitudinal momentum matrix elements taken between the conduction and valence 
band edge Bloch functions^, and mf and mf are the remote-band contribution to the 
transverse and longitudinal band edge effective masses, respectively. For electrons and holes, 
these band edge effective masses can be expressed as mf t = [l/mf t + 2Pj 2 t /m 1 E g }~ 1 and 
m\ t = [l/m"j r t + 2Pj 1 t /m l E g }~ 1 , respectively. In each valley, the z axis in Eq. Q is directed 
toward the L-point of the Brillouin zone, e.g. along the (111) direction of the cubic lattice. 
As a result, for each of the four valleys, the z axis will point in different directions. 

Although the Hamiltonian of Eq. ([3]) has cylindrical symmetry with respect to, e.g., 
the (111) crystallographic direction, this direction may not coincide with the NR growth 
direction. For a description of NR electronic and optical properties it is convenient to 
use coordinates connected with the latter direction instead, even though the cylindrical 
symmetry of the Hamiltonian is generally broken. In PbS and PbSe, the small anisotropy 
of conduction and valence bands allows us to treat deviations from cylindrical symmetry 



perturbatively. The Hamiltonian can be written H = H + H an , where the cylindrically- 
symmetric part Hq is 




The modified band parameters are 



m 




2 9 



sin 2 9 



P z = P t sin 2 9 + Pi cos 2 9 




-\ sin 2 9 + cos 2 9 (6) 



raj 



(5) 



where 9 is the angle between the growth axis and the (111) direction. The anisotropic 
part of the Hamiltonian is given in Appendix |A] Note that Eq. (J4]) has a form identical to 
Eq. ([3]), but the z axis is now directed along the growth axis. For arbitrary orientation of 
the growth direction, there will be four different angles 9 for each of the four valleys, and 
therefore four different sets of modified band parameters defined in Eq. ([5]). As a result, 
each valley will have unique electronic structure. 

The energy spectra associated with the different valleys become degenerate when the 
growth direction leads to identical values of 9 for them. The highest degree of degeneracy 
is reached when the growth direction is along the (100) crystal axis. In this case all four 
valleys have the same 9 ; cos 2 9 = 1/3 , which results in P± = P z and mi = m z in Eq. (T4j) . 
All of the spectra are degenerate. 

The anisotropic part H an of the full Hamiltonian can be considered as a perturbation 
if \Pi — P t \ <C Pi + Pt and \l/mf — l/mf | -C l/fnf + l/mf . The first-order corrections 
to the solutions of Hq caused by H an vanish in the 2-fold Kramers-degenerate subspace at 
each energy level. As a result, only second-order perturbation theory gives corrections to 
the energy levels. We will neglect these corrections from this point on, although an example 
higher-order calculation appears in Appendix |A] 

III. ENERGY SPECTRA OF ELECTRONS AND HOLES IN PBSE NANOWIRES 

The first step in our modeling process is to find the energy spectra of ID subbands of 
infinitely-long cylindrical nanowires, temporarily ignoring the Coulomb interaction. The 
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cylindrical symmetry of the Hamiltonian of Eq. (J4j) allows the solutions to take the form 

( 7^0)exp(z(n- 1/2)0) \ 
iWZ(p) exp{i{n + 1/2)0) 
7^(p)exp(i(n- 1/2)0) 
Vft2(p)exp(z(n + 1/2)0)/ 



exp(i/c^z) 



(7) 



where is the azimuthal angle, n = ±1/2, ±3/2, ±5/2, ... is the total angular momentum 
projection on the nanowire axes defined by the operator J z = —id/d<f> + S z , hk z is the 
momentum along the nanowire z axis, and p = \J x 2 + y 2 is the radial coordinate in the 
plane perpendicular to the NW axis. The chosen phase of each component of the function 
W l {k z ) allows the radial functions 7£™(p) to be pure real. Substitution of Eq. ([7]) into Eq. 
(00) yields the system of differential equations that defines these functions: 



Ir 



2m_! 

h 2 

2m", 



-A 



n-l/2 



m m 



fth P 

-K(p) 



m 

hP 



hP± 

m 



-D 



n+l/2 



h 2 



m 



m 



2m 

_^ 

2m 



-A n _i/ 2 

A n+ i/2 



K(p) 









(8) 



where a± = E g /2 ± E + h k 2 / (2mf) . The differential operators 



D 



± 



=F 



<9p 



m 
P 



(9) 



are the raising and lowering operators D^J m (kp) = kJ m ±\(kp) for the Bessel functions 
Jm(kp) with integer index, and the operator A m = D m+l D+ n = —(l/p)(d/dp)p(d/dp) + 
m 2 /p 2 . 

It is easy to check using the raising and lowering properties of the operators that 
the radial eigenfunctions of Eqs. ([ED should take the form 



K n 2 (p) 
T^(P) 



(c 

. Jn—1/2 
C2Jn+l/2{k p p) 

CzJ n -i/ 2 {k p p) 
\C 4 J n+1/2 {k p p) J 



(10) 
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Substitution of this into Eqs. dHJ yields a 4x4 system of linear equations for the coefficients 
Ci.2,3,4 • Setting the determinant of this system to zero produces the relation between the 
quasi-momentum k p and the energy of electrons or holes E : 



h 2 k 2 p = -a(E) ± vM^) 2 + P(E) 



(11) 



where 
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2m+ 2 / \ 2m^ 2 y m^ 
From Eq. (TTTj) it is clear that fc^ can be positive or negative. The negative value of k 2 results 



in an imaginary k p = i\ p , with X p defined by Eq. fflTl) as h 2 X 2 = a(E) + a/ 'a(E) 2 + ^(-E) . 
The complex arguments in Eq. (TTTJ]) are then simplified by replacing the Bessel functions 
Jm{iX p p) with the modified Bessel functions I m (X p p) using the relationship J m {iX p p) = 
i m I m (X p p) . For each value of k 2 , there are two independent solutions of the 4x4 linear 
system for the coefficients 6*1,2,3,4 . These two solutions can be chosen such that either 
C3 = or C4 = , which allows the remaining coefficients Cj to be found. Taking into 
account the positive and negative value of k 2 , there are four independent solutions for each 
energy E . 

The energy spectrum is determined by the boundary conditions at the NW surface. The 
boundary conditions are defined on all four components of the wave function, which inside 
of the NW can be always written as a linear combination of the four degenerate solutions 
discussed above 
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where 

and a ,b , c , and are determined by the boundary conditions. 

For NWs with an impenetrable surface, the standard boundary conditions require each 
component of the wave function defined in Eq. (TT3]) to vanish, leading to TZ^R, k z ) = , 
where i = 1,2,3,4 and R is the NW radius. These four equations define the 4x4 system 
for the a, b, c, d coefficients. Requiring the determinant of this system to be zero yields the 
following implicit equation for the ID energy bands for angular momentum n , and as a 
function of the parameter k z : 

t,2p2(-p _ ri \2 , p2(h,2-p2 _ \2p2\ 

K\ [(il)V-) 2 - (i n S(Jl) 2 ] + -^-^ — p 2 rr x p k) i n -i n + J-Jl = o , (is) 

i t J- k± A 

where we use the notation J± = J n ±i/2{k p R) and I± = I n ±i/2(^ P R) ■ 

After determining the energy from Eq. (fT5l) . the wavefunctions can be constructed from 
Eq. (|T3|) . with only the normalization undetermined. We will use the following notation 
for normalized eigenf unctions: \l/g' fc and for the electron and hole levels given by Eq. 
(TL5l) . correspondingly, where k = 1,2,3... is the index of the ID subband with angular 
momentum n , and \ ty™' k \ 2 p e dp e 2ir = jj^ \ ^' k \ 2 p h dp h 27r = 1 . 

Using Eq. f|T5|) we calculated the energy levels for a 4-nm PbSe NW with various growth 
directions. The energy band parameters of PbSe which we used in this calculation will 
be described in a later section. The effective energy gap of the NW, which is the energy 
distance between the top of the highest ID sub-band of the valence band and the bottom of 
the lowest ID sub-band of the conduction band, impacts many material properties. Figured] 
shows the effective energy gap for all four valleys as a function of the growth direction of the 
nanowire. Because the plot is calculated along high-symmetry directions in the Brillouin 
zone, the degeneracy of the four valleys is never completely lifted. Without any intervalley 
coupling, each of these energy gaps would have separate optical absorption and emission 
peaks associated with it. 

Figures [2^ and [2b show the dispersion of the several lowest ID subbands of the conduction 
and valence bands in NWs grown along the (111) and (100) directions, respectively. NWs 
grown along (111) have one valley oriented parallel to the growth direction and the other 
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<111> <100> <110> <111> 

Rod Growth Direction 



FIG. 1. Energy gaps of a 4-nm diameter PbSe NW at each of the four valleys as a function of 
the growth direction of the NW (red lines). The numbers indicate the valley degeneracy of the 
energy gaps. Dashed grey lines are the same energy gaps after accounting for the self-Coulomb 
interaction, described later in the text. 

three valleys oriented at the equal angles 9 = 71° from it. For the (100) NW, all four 
valleys are at the same angle 9 = 55° from the growth direction. It is clear that both 
the band-edge energies and the effective masses of the ID subbands depend strongly on the 
growth direction. 



IV. DIELECTRIC CONFINEMENT 



The optical properties of all semiconductor nanostructures are controlled by the strength 
of the Coulomb interaction between the electron-hole pair participating in the emission 
and absorption of photons^. Compared to the screened Coulomb interaction in a bulk 
crystal, the interaction is usually enhanced because the electric field of the electron and hole 
localized inside the nanostructure penetrates into the surrounding medium, which commonly 
has a dielectric constant smaller than that of the semiconductor. In addition, any charge 
in the vicinity of this interface polarizes it. In the case of a flat interface, for example, this 
polarization can be described easily using an image charge that interacts with the primary 
chare; In the case of small external dielectric constant the interaction is repulsive. This 
repulsive potential in nanostructures of any shape leads to an additional confinement of 
carriers, which is referred to as dielectric confinement. 



9 



-2-10 1 2 

k z [nm"'] 



-2-1012 
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FIG. 2. ID band structure of a 4-nm PbSe NW for the cases of the axis along the directions (a) 
(111) and (b) (100) . The bands are labeled by the angle between the considered valley and the 
rod growth axis and also by their multiple valley degeneracy, up to a maximum of (x4). In (b), the 
individual subbands are labeled using notation adopted from molecular physics: kX^ for the k th 
electron or hole level of certain symmetry with total z angular momentum n , where X = £ , II , 
A ,. . . , is used for \m\ = , 1, 2,. . . , respectively, where m is the angular momentum projection 
of the conduction (valence) band component of the wavefunction of the electron, ' e ', (hole, ' h ') 
state. In (a), the order of the levels is the same, and the labeling is suppressed for clarity 

To model these effects in NRs and NWs, the analytic potential for two charges in an 
infinite dielectric cylinder U(r e ,rh)— is used. It was shown previously-^ that this approx- 
imation works well as long as the rod length is larger than the size of the exciton. The 
potential naturally divides into four terms^i: the unscreened direct interaction of the two 
charges Ud , the modification of this interaction due to the image effects of the solvent U s , 
and the two self-interactions of each charge with its own image U e and Uy t : 



K s K m (Ru) (J m _i(Hit) + I m+1 (Ru)) + K m I m (Ru) (Km^Ru) + K m+1 (Ru)) 
and where k s and the optical dielectric constants of the bulk semiconductor and 

the surrounding medium, respectively. I m and K m are the modified Bessel functions of the 



U(r e ,r h ) 



e 2 /(K s \r e -r h \) - eV s (r e ,r h ) + eV s (r e ,r e )/2 + eV s (r h ,r h )/2 
U d (\r e -r h \) + U s (r e ,r h ) + U e {r e ) + U h (r h ) 



(16) 




(17) 
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first and second kind. For PbSe we will use k s = 23 , and for the medium, if not explicitly 
stated otherwise, n m = 2 throughout this work. 

The self-interaction terms U e (r e ) and Uh{rh) always contribute to the energy of each 
electron and hole subband calculated in Section [TTTJ In narrow NWs and NRs, where the self- 
interaction energy is smaller than the confined energies, this contribution can be calculated 
perturbatively for electron and hole levels, respectively: 

Kie = f Pedped^m^U^) , = J Ph d Ph d^l ' \ 2 U h (r h ) . (18) 

The self-interaction terms E^ {e and E™ c ^ h increase the energy of all electron and hole ID 
subbands and consequently the effective energy gap in nanowires. The perturbed electron 
and hole subbands with n — n' — 1/2 and k — k' — 1 are shown in Fig. [TJ 

In addition, in narrow NWs and NRs one can used an adiabatic approximation of the 
Coulomb interaction^ 2 -^, which replaces the three-dimensional potential of electrons and 
holes of Eq. f|T6|) by a one-dimensional Coulomb potential that describes their interaction 
along the NW/NR axis. The adiabatic potential is obtained by averaging the potential over 
wave functions fy™' k and )k of the corresponding electron and hole subband. Averaging 
the first two terms of Eq. (fT6l) results in the ID adiabatic potential 

V^\\z e - z h \) = j Pe dp e d<f) e J p h dp h dif> h \^ k \ 2 \^' k '\ 2 {U d (\r e - r h \) + U s (r e ,r h )) , (19) 

which describes the interaction of electrons and holes occupying different subbands. This 
adiabatic potential is a function of the electron and hole separation, \z e — Zh \ , only. One can 
show that at large distances \z e — Zh\ ^> R it takes the form of a one-dimensional Coulomb 
potential with the dielectric constant of the surrounding medium, V™£ ~ — e 2 / (K m \z e — Zh\) ■ 
The adiabatic potential for the ground electron and hole subbands with n — n' — 1/2 and 
k — k' — 1 is shown in Figj3j 

V. ID EXCITONS IN PBSE NANOWIRES AND NANORODS 

The attractive ID potential described by Eq. f|T9l) creates a series of one-dimensional 
exciton states for each pair of electron and hole subbands (n, k) and (n' , k') . The effective 
masses of electrons and holes along the NW axis m™ ,fc and m™' ,fe ' at the bottom and the 
top of each subband, correspondingly, is determined by Eq. (TT5I) . This allows us to write a 
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FIG. 3. Points show the effective binding potential, V^L ^ , between an electron and a hole 
occupying the ground one dimensional subband n = n' = 1/2 and k = k' = 1 as a function of 
their separation, calculated for a 4-nm radius PbSe NW. The solid line shows the approximation 
of this dependence by the Elliott & Loudon effective potential described by Eq. (|2ip 

one-dimensional Schrodinger equation for these ID excitons: 



h 2 d 2 



h 2 d 2 

2M^' dZ 2 



*iD+£c:f(2)*iD=<r * 1D 



2fJi n'k' g z 2 ' ( 20 ) 

where we introduce the electron-hole separation, z = z e — z^ and the exciton center-of-mass 
coordinate Z = [m n /z e + m^' k ' z h )/(m2' k + m^' k ') . = m n /m n ^ k ' j{m n e > k + m"'' fc ') is the 

reduced mass and M™'^ = m™' k + ri\ ,k is the total effective mass of the ID exciton. Impor- 
tantly, the exciton binding energy in this equation is calculated relative to the distance 
between the bottom of the (n, k) conduction subband and the top of the (n', k') valence 
subband, assuming the self-interaction energy terms E^ k ie and E r s l c { k h are already taken 
into account. The solution of Eq. fl20|) can be separated into \[> 1D (z, Z) = tpi D (z)^ cm (Z) . 
The wave function ^id(^) describes relative electron- hole motion and gives the spectrum of 
ID excitons. The second component, ^ cm (Z) , describes the exciton center of mass motion, 
and in the case of an infinite NW \l/ cm (zT) ~ exp(iKZ) , where hK is the exciton momen- 
tum along the NW axis. This replaces the second term in Eq. ([20]) by the exciton kinetic 
energy, h 2 K 2 /2M^'. 

Equation ( 120]) allows us to numerically calculate the energy spectrum of ID excitons 
created from any pair of electron and hole subbands. In this paper, we will be interested 
primarily in the spectrum that arises from the lowest electron and hole subbands 1S^ 2 and 
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l£jy 2 , and we will use the approach suggested by Elliott & Loudon^ to describe the spec- 
trum of one- dimensional excitons in a strong magnetic field. They suggest approximation 
of the one-dimensional adiabatic potential by an effective one-dimensional potential, which 
has well-known Schrodinger equation solutions, 



U eS (z) 



Ap eS e 2 



(21) 



«ro(|«| + Pes) Km(\z\ + Pcff) 2 ' 

where p C Q and A are the two fitting parameters. The medium dielectric constant K m is 
used in Eq. ( 12TT) so that the correct asymptotic form of the potential is maintained. For a 
4-nm PbSe NW immersed in a medium with n m = 2 , the numerically-calculated effective 
potential is described very well by the potential U e g with p e s = 5A9R and A = 2.73 , as 
seen in Fig. |3j The slight dependence of these fit parameters on NW size is shown in Fig. 
and the much stronger dependence on K m is shown in Fig. Hb. 
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FIG. 4. Fitting parameters used in the effective potential described by Eq. (|2ip in PbSe NWs of 
various radius and medium dielectric constant. The parameter is plotted vs. (a) nanowire diameter 
with K m = 2 (b) medium dielectric constant with R = 2 nm. 

The energy spectrum and eigenfunctions of Eq. (j20l) with effective attractive potential 
U e s(z) can be obtained analytically. The eigenfunctions of each ID exciton level, i/} a {z) , 
can be written as^ 1 ^ 



ip a {z > 0) = a 1 W a + p) + a 2 M a ^i VT ^ A ^{z + p) 

i) a {z<V>) = ±iP a {\z\) 



(22) 
(23) 



where W a ,p(x) and M a ^(x) are the Whittaker functions, z = 2z/(aoa) , p = 2p c s/(a a) , 
a = h 2 K, m /(p\^'\e 2 ) is the effective Bohr radius of a ID exciton, and a\ and a 2 are 
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arbitrary coefficients. The sign of Eq. (123|) is " + " for an even eigenfunction and " — " for 
an odd one. The coefficients a\ , 0,2, and parameter a in Eq. ( 122]) as well as the exciton 
binding energy: 



h 2 



are determined by the boundary conditions. 

There are two boundary conditions to impose on the solution in Eq. f[2"2"l) : one at z = 
z e — Zh = ±L and one at z = . We first consider infinite nanowires; the effects of finite 
length will be treated in the following section. In this case, the first boundary condition is 
satisfied by letting a 2 = , because M a _i^/ 1 _ 4Aa - (\z\+p) diverges as \z\ — > 00 . The second 
boundary condition, requiring ip a (z) to be either an even or odd function of z , determines 



(24) 



a and the energy spectrum of the exciton. It was shown in Refs. |32| and[33[ that for excited 
doubly-degenerate exciton states, a takes almost-exactly integer values a = 1,2,3,... and 
hat a — > for ground states with decreasing exciton transverse radius. Following Refs. 



32 



and 



33( we use Eq for the ground exciton binding energy. 
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FIG. 5. Coulomb energies calculated for (a) K m = 2 with varying R and (b) R = 2 nm with 
varying K m ■ Lines are the sum of the electron E^'^ and hole ^ self interaction energies 
(red circles); the electron- hole binding energy eo (blue triangles); and their total (black squares). 

Figure shows the calculated binding energy of the ground exciton state £0 an d the 
Coulomb self-interaction energies E^'l and '\ of electrons and holes from the ground 
ID subbands lS^ ■ The binding energy decreases dramatically with NW radius or external 
dielectric constant. The exciton binding energy in the narrowest NW surrounded with 
K m ~ 2 — 3 reaches values > 300 meV. 
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Surprisingly, however, the binding energy is almost exactly compensated by the electron 
and hole self-interaction terms, which leads to practical cancelation of most effects connected 
with the small dielectric constant of the surrounding medium. Because of this cancelation, 
the optical transitions between ID subbands will be determined primarily by the energies 
calculated in Section IHII This result has important practical consequences. For example, 
the linear optical spectra of PbSe NWs will not be sensitive to the dielectric constant of the 
surrounding medium. 

This cancelation is well-known in spherical semiconductor NCs. The exact cancelation 
of these three terms was shown for parabolic valence and conduction bands in Ref. 34]. 
This is because in a parabolic-band approximation the wave function of electrons and holes 
are identical and depend only on the NC radius. As a result the electron and hole charge 
distributions exactly compensate each other at each point in the NC. If there is no local 
charge in the NC, there is no electric field outside of the NC, and the external medium does 
not affect the optical properties. This cancelation is nearly exact even when the electron 
and hole masses are different^. 

The cancelation of the Coulomb energies in the ground exciton of PbSE NWs can be 
attributed to a similar charge compensation. The mirror symmetry of the conduction and 
valence bands in PbSe makes the wave functions of the electron and hole transverse motion 
nearly identical. The similar values of effective masses along the NW axes also makes the 
electron and hole contributions to the ID exciton wave function identical. It is interesting to 
note here that because of the large binding energy, the electron and hole in the exciton are 
remarkably tightly bound, with average separation only slightly larger than the NW radius. 



Fig. [6] shows the average separation, calculated as \J ((z — z) 2 ) , as a function of radius, 
with inset showing the wavefunction ifiiz, for the case of R = 2 nm. One can see that the 
average electron-hole separation in the exciton is an order of magnitude smaller than the 
46 nm Bohr radius in bulk PbSe. Further calculations show that this unusual increase in 
the strength of the binding is due entirely to the ID shape of the NR, and is only weakly 
affected by the dielectric contrast. For the weakest dielectric contrast when K m = k s = 23 , 
the average separation increases slightly to ~ 4 nm, still much closer to the 4-nm diameter 
than to the Bohr radius. 
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FIG. 6. Dependence on PbSe NW radius of the average (r.m.s.) separation of the electron and 
hole in the ID exciton. Inset shows the square of the ground exciton wavefunction IV'idI 2 for a 
NW with 2 nm radius. 

A. Finite length effects 

For a nanorod, which has finite length, the relative and center-of-mass (CM) motions 
of the electron and hole can never be completely separated. If the NR is much longer 
than the radius of the ID exciton, one can still approximately separate variables to create 
effective boundary conditions for the exciton CM motion. No other boundary condition 
(BC) is needed for the exciton separation coordinate, because the assumption is that the 
tightly-bound wavefunction is already zero well before any additional confinement is felt. 
On the other hand, the CM motion can be considered as the motion of a free particle 
confined in a ID box of length L . If the box is much larger than the exciton radius one 
can apply the standard boundary conditions on \I/ cm to obtain the well-known spectrum 
E cm (l) = h 2 7T 2 P 1 1 (2M\^\L 2 ) , where / is the level number. 

Even though this CM boundary condition makes intuitive sense, it is difficult to justify, 
because the true BCs are for the electron and hole individually. To test our assumption, 
we calculated the CM energies numerically by solving the two-particle Schrodinger equation 
with the correct impenetrable boundary conditions on the electron and hole individually. 
Details of the calculation are in Appendix The numerically calculated wavefunctions 
and energies were best matched to those obtained for a free particle with an effective mass 
of the exciton which is confined in the ID box of length L cm = L — R. The existence 
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of such a simple expression is probably connected with the approximately-equal effective 
masses of the electrons and holes and their small separation in PbSe NRs. The first few 
numerically-calculated energy levels are shown in Fig. [7J along with the analytic energies 
E cm = H 2 l 2 7T 2 / (2My^L 2 m ) for various confinement lengths L cm . This modified CM length 
works well for all rod sizes studied, as long as the NR aspect ratio is > 2 . 
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FIG. 7. Numerically-calculated energies for the lowest few exciton states in a 4 x 20 nm PbSe NR 
(black circles). The lines are the energies from the analytic model using two different confinement 
lengths for the center of mass. 
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B. Oscillator strength of the interband optical transitions 



The decrease of the electron-hole separation within a ID exciton leads to a dramatic 
increase of the optical transition strength. It was shown by Elliott & Loudon^ that the 
oscillator strength of practically the entire spectrum of ID excitons becomes concentrated 
in the ground exciton state. The expression for the transition strength in PbSe NRs can 
be obtained by combining the results derived for PbSe NCs^ and CdSe NRs^. The total 
oscillator strength to tai can ^ e wr hten as a product O to tai = 0±0\\ , where the tranverse 
oscillator strength is^ 
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with hio the total energy of the optical transition. We have neglected the second term 



from Ref. 



24j . as it is negligible except for very small NRs, where the envelope function 



approximation likely breaks down anyway. The oscillator strength of the ID exciton^. is 



On = |^ 1D (z = 0)p 



dZV cm (Z) 



(26) 



where we normalize the ID exciton wave function such that J^ L dz J Q L dZ\i/ji D (z)^/ cm (Z)\ 2 = 
1 . 

The transverse oscillator strength provides the selection rule that there is no change in the 
z-component of the angular momentum, An = , while the longitudinal component focuses 
the oscillator strength into the ground exciton state. This is because optical transitions 
are only allowed to the even states of the exciton CM motion with I = 1,3,5... , and the 
oscillator strength decreases as l/l 2 . Even the second allowed transition will be 9 times 
weaker than the lowest transition. This has practical implications for the optical absorption 
spectra. Even though the density of allowed transitions increases dramatically with energy 
in NRs, most of the oscillator strength is concentrated in the lowest-energy transition for 
each pair of NR subbands. Thus, isolated peaks should still be observable in experimental 
spectra. 

VI. EXPERIMENT 

A. Synthesis and characterization of colloidal PbSe nanorods 

Although the synthesis of lead salt nanowires was reported several years a,gcM^-, the 
fabrication of high quality lead-salt nanorods with small diameter has proved challenging. 
PbSe NRs were synthesized with noble metals as seeds^, but the resulting NRs did not have 
good optical spectra. Some high-quality NRs have been reported, but the syntheses were 
too challenging for us to reproduce^— . A simple synthesis for high-optical-quality PbSe 
NRs was recently demonstrated^, and the properties of these NRs will be compared to the 
theoretical results. 



Following Ref. [421 ]. the NR synthesis was carried out using standard Schlenk-line tech- 
niques under dry nitrogen. Tris(diethylamino)phosphine (TDP, Aldrich, 97%), oleic acid 
(OA, Aldrich, 90%), 1-octadecene (ODE, Aldrich, 90%), squalane (Aldrich, 99%), amor- 
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phous selenium shots (Se, Aldrich, 99.999%), and lead(II) oxide (PbO, Aldrich, 99.9%) were 
used as purchased without further purification. Anhydrous ethanol, chloroform, acetone, 
hexane, and tetrachloroethylene (TCE) were purchased from various sources. To prepare 
1.0 M stock solutions of TDPSe, 7.86 g of Se was dissolved in 100 mL of TDP. 

Typically, 0.22 g of PbO was dissolved in 5 mL of squalane in the presence of 1 mL 
OA. (Squalane can be replaced by ODE.) After drying under nitrogen at 150 C for 30 min, 
the solution was heated to 170 C and 3 mL of a 1 M TDPSe solution in TDP was injected 
under vigorous stirring. Once the reaction finished, the reaction mixture was cooled to room 
temperature using a water bath. The crude solution was mixed with hexane and precipitated 
by ethanol. The precipitated NRs were isolated by centrifugation (at 5000 rpm for 3 min) 
and redispersed in chloroform or other organic solvents. Size-selective precipitation can be 
carried out to obtain better monodispersity of NRs samples using chloroform/acetone or 
other solvent /nonsolvent pairs. 

The size of the synthesized NRs was determined from transmission electron microscopy. 
In-plane powder X-ray diffraction shows that the NRs grow along the (100) direction^. 
Absorption was measured on a Shimadzu UV-3101PC spectrophotometer at room temper- 
ature. Emission spectra were recorded at room temperature with an infrared fluorimeter 
equipped with a 200-mm focal length monochromator, a single mode fiber coupled laser 
source (S1FC635PM, 635 nm, Thorlabs, Inc) as the excitation source, and an InGaAs 
photodiode (New Focus Femtowatt model 2153). Fluorescence lifetime was measured us- 
ing an InP/InGaAs PMT (Hamamatsu H10330A-75) with 120-fs excitation pulses from a 
Ti:sapphire regenerative amplifier (Spectra-Physics Hurricane) with 1 kHz repetition rate. 
NRs were dissolved in tetrachloroethylene (TCE) for all measurements to avoid spurious ab- 
sorbance in the near-IR. Quantum yield measurements were performed using an integrating 
sphere, with the method described in Ref. 43]. 



B. Absorption Spectra 

First, we will highlight the qualitative differences between the absorption spectra of NRs 
and spherical NCs. Figure El shows the absorption spectrum of 3.3 nm diameter x 12 nm 
length PbSe NRs along with that of 4.4 nm diameter spherical NCs, chosen to have a nearly 
identical first absorption peak. The spectrum of the NRs has fewer obvious features than 
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the NC spectrum. The first peak in the NR spectrum has a broad high energy side, even 
though its narrower low energy side is nearly identical to that of the NCs (inset of Fig. [8]). 
Both of these observations indicate the presence of more densely-spaced transitions in the 
NR spectrum, which have the effect of smoothing out the peaks. Interestingly, the second 
NC peak appears where there is a dip in the NR spectrum. 
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FIG. 8. (a) Absorption spectra of PbSe NRs (black line, vertically offset for clarity) and spherical 
PbSe NCs (red line) are compared. The inset shows detail of the first peak, (b) Emission spectra 
and fluorescence decays measured at the emission peak (inset) of the same two samples. 



The broadening of the NR absorption peak seen in Fig. [8] is connected with the dispersion 
of NR diameter and length. Our best PbSe NR samples have around 5% size distribution 
in radius, but a much larger 20% in length. This large length polydispersity will blur out 
many of the NR transitions in an ensemble, except for those that are roughly independent of 
length — specifically, the lowest energy exciton for each pair of NW subbands. Fortunately, 
this is also the transition predicted to have the largest oscillator strength. As we have 
shown above, the energies of the optical transitions of the ground exciton states practically 
coincide with the energies between non-interacting electron and hole subbands, even though 
their respective wave functions differ greatly. This greatly simplifies the interpretation of 
the absorption spectra of NRs. 

We performed second-derivative analysis on the absorption spectra to determine the 
transition energies accurately. To avoid the problems inherent in this method^, only the 
peaks in the second- derivative spectra that correspond to obviously-visible peaks in the 
measured spectra were used. NRs produced by our first syntheses showed instability in 
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solution and would slightly aggregate during the absorption measurement. This adds a 
moderate scattering background, so only the absorption peak location is recorded for these 
samples. NRs synthesized more recently are more stable, and at least four peaks can be 
discerned, with an additional peak in the three samples with narrowest size distribution. 
Fig. [Hh. has an example measured spectrum of a 3.9 nm diameter x 17 nm length PbSe 
NR that shows all five peaks, and the locations of all measurable peaks from all samples 
are shown in Fig . [9b. The measured peaks are plotted vs. D~ 3 / 2 following the similar 
graph in Ref. 45j. This power of the diameter is chosen to make the trend linear over the 
measured range, allowing rough extrapolation to bulk as D~ 3 ^ 2 — > . In this manner, the 
peaks originating from the L-point and £ -point are easily distinguished. 
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FIG. 9. (a) Example absorption spectra of a 3.9 x 17 nm PbSe NR. Inset shows the same data, 
but on a scale where the 5th peak is visible, (b) Peaks in 2nd derivative spectrum as a function of 
NR diameter (symbols), calculated allowed transitions (grey lines), simple parabolic effective mass 
calculation around the £ -point (dashed grey line), and linear fits (colored dashed lines.) 

Quantitative theoretical description of the size-dependent absorption spectra of PbSe 
NRs shown in Fig. [9] requires a set of 6 room temperature energy band parameters for 
this semiconductor: mf , mf , and . The parameters extracted from low temperature 
cyclotron resonance and interband magnetooptical experiments in bulk PbSe^ describe quite 
well the average two-dimensional effective mass of electrons and holes at the bottom of the 
conduction band and the top of the valence band, respectively. The fitting procedure that 
gives this set is not sensitive, however, to the separation of and the 2Pf t jm\E g 

terms, and describes well only the sum of these terms, because the all measurements are 
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conducted a the narrow energy range comparable with the PbSe energy gap. This procedure 
is also not very sensitive to the anisotropy of the carrier energy spectra, because a magnetic 
field averages out the 2D motion of electrons and holes. On the other hand, in order to 
predict nanocrystal energy levels quantitatively, both the separation of components of the 
effective masses and the band anisotropy are crucial. Finally, the energy band parameters 
are expected to be temperature dependent. Thus, we conclude that parameters inferred 
from cyclotron resonance and magneto-optical measurements might not describe the energy 
spectra of NRs and NCs measured at room temperature. 

In principle, spatial confinement of carriers in nanostructures provides a more-sensitive 
way to determine the energy band parameters, due to the large modification of the energy 
spectra of confined carriers. With this motivation, we used the previously-measured ab- 
sorption spectra of PbSe NCs in Refs. 
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47H51| and extracted room-temperature band 



parameters using a global fitting procedure. Importantly, this new set of parameters not 
only quantitatively describes the low-energy transitions of PbSe NCs, but may also help 
resolve the long-standing controversy over the symmetry of the second peak in the NC ab- 
sorption spectra (see Appendix IDl). These band parameters (Table IVI Bj) are used in all 
graphs presented in this work. 



Name 


Ref. [46] 


Best Fit 


Aniso. ratio 


Rcf. [46] 


Best Fit 


mf/m 


0.29 


0.59 


mt /m^ 


1.28 


1.6 


m// /mo 


0.27 


0.79 


mT /m^ 


3.53 


1.6 


2P t 2 /m 


3.6 (eV) 


4.25 (eV) 




1.82 


3.0 



TABLE I. Energy band parameters that provide the best fits to the room temperature data from 
PbSe NCs. The left columns show the transverse band components, while the right columns show 
the ratio of transverse to longitudinal components. 



The theoretical size dependence of the optical transitions in PbSe NRs is calculated within 
our 4 band model and shown in Fig. [9b by solid lines. The lowest two transitions agree 
well with the theory. The third predicted transition is not observed, possibly owing to its 
proximity to other strong transitions in our NR samples. The third and fourth peaks are 
strong transitions that do not appear to be associated with the L-point. Their energies 
extrapolate back to the S -point energy. The third peak is fit well by the same parabolic 
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band model used to model spherical PbSe NCs, and thus we assign this transition to the 
lowest-energy excitonic state at the E point. This line was calculated for both spheres and 
rods with m| = m\ = 0.45m and E g (E) = 1.65 eV. Without more-detailed knowledge of 
the band structure there, we cannot predict the excited states with any accuracy. Thus, the 
identity of the fourth transition cannot be determined, but as the energies approach the same 
1.65-eV bulk value, it is reasonable to tentatively attribute it to a higher-energy exciton from 
the E point. Finally, the fifth peak was perhaps the strongest in the absorption spectra, but 
showed no size dependence. We tentatively ascribe this to a metal-complex transition on the 
surface of the nanocrystal based on its proximity to absorption peaks of Pb(II) complexes^. 
The identities of these transitions are summarized in Table IVIBl 



Label 


Assigned Transitions 


Pi 


1S l/2 1S l/2 


P 2 


m 3/2 m 3/2 ^d m h 1/2 m? /2 


Ps 


S-point ground state 


Pa 


S-point excited state (?) 


P 5 


Surface metal complex mode 



TABLE II. Transitions observed in the absorption spectra of PbSe NRs. 

The fluoresence spectra and decays (Fig. [8b) are nearly identical for NCs and NRs, with 
a slightly larger Stokes shift in the NRs along with a slightly broader peak. The ensemble 
quantum yield of the nanorods is around 15%, around half that of the NCs. This might 
indicate that the radiate lifetime of the rods is longer than the that of the NCs, but it is also 
possible that the QY reflects an ensemble with 15% emitting and 85% non-emitting rods. 

Two effects would be expected to modify the radiative lifetime in nanorods. First, be- 
cause the radiative lifetime is inversely proportional to the oscillator strength, the increased 
electron-hole correlation in NRs should decrease the lifetime compared to NCs. Second, the 
effect of screening is reduced in NRs, which is believed to be the cause of the long lifetime 
in PbSe NCs^. Approximating the NR as a dielectric prolate spheroid, the screening will 
substantially decrease along the rod axis, while slightly increasing along the other two axes, 
with an overall effect of a reduction in screening of the lifetime. Compared to a spherical 
NC of the same diameter, the larger oscillator strength and the reduced screening should 
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each produce about a factor of 3 reduction in lifetime in NRs with typical aspect ratios. To- 
gether this amounts to almost an order of magnitude reduction, and should be measurable 
even considering other sample-related uncertainties. However, the measured lifetime (Fig. 
[Sb) is nearly identical in NCs and NRs. This discrepancy is not understood. It might be 
explained by a dark ground exciton state that controls the photoluminescence decay in PbSe 
NRs and NCs, with the same activation mechanism in both structures. To be thorough, 
the nonradiative rate must be determined, and completing this along with exploring this 
phenomenon is a topic of future work. 



VII. DISCUSSION AND CONCLUSIONS 



Our model of the electronic structure of lead-salt NRs is based on the 4 band k ■ p 
Hamiltonian suggested in Ref. 0, using the standard boundary condition of a vanishing 
envelope wave function at the NR surface. All calculations are conducted within a cylindrical 
approximation. To use this model for description of various properties of NRs or NWs, one 
needs to know a set of the 6 temperature-dependent band parameters that describe a specific 
bulk lead-salt semiconductor. For the PbSe NRs studied in this paper, we extracted the 
set of room-temperature parameters from analysis of the size- dependence of previously- 
measured room temperature absorption spectra of spherical PbSe NCs. 

The most significant conclusion of this work is that the fundamental excitations in PbSe 
NRs are one- dimensional excitons under each pair of optically coupled electron-hole sub- 
bands. The binding energy of the ground exciton state, which accumulates the most oscilla- 
tor strength, increases with decreasing NR thickness and reaches 400 meV in the narrowest 
rods. Surprisingly, the large binding energy of the exciton is almost exactly compensated by 
the self-interaction of electrons and holes with their own images, which makes the energies 
of the optical transitions nearly independent of the solvent dielectric constant. Although the 
finite length of NRs affects the spacing between excited exciton states, it has a negligible 
effect on the energy of the exciton ground states. 

With the set of PbSe band parameters extracted from spherical NC absorption spectra 
(Table [VI Bp . the model presented here describes the absorption spectra of PbSe NRs, and 
potentially resolves some troublesome aspects of k ■ p theory of spherical PbSe NCs. The 
energy of the optical transitions to the exciton ground states calculated within a cylindrical 
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approximation match the two lowest-energy transitions observed experimentally. Although 
the effect of anisotropy in important for description of the absorption in spherical PbSe NCs, 
it is diminished in NRs (see Appendix |A]&[D]), and the energy of the first two transitions is 
unaffected by it. 

The absorption spectra of PbSe NRs have another remarkable feature. The size depen- 
dence of the third and fourth absorption peaks is strong evidence that they originate from 
the £ point of the Brillouin zone. Similar states connected with the £ point were ob- 
served previously in the absorption spectra^ and in the hot carrier dynamics^ of spherical 
PbSe NCs. These observations provide clear experimental evidence that even in the smallest 
nanostructures, wave functions from distinct critical points ( L and £ , in this particular 
case) are not mixed if both their corresponding conduction band minima and valence band 
maxima are well-separated energetically. This provides strong justification for the applica- 
bility of our multiband effective mass approximation in such small nanostructures. A large 
energetic separation of L and £ band edges is supported theoretically by recent ab ini- 
tio calculations^, which for PbSe predict larger than 500 meV energy separation for these 
extrema, in both the valence and conduction bands, in contradiction with similar earlier 
calculations, which placed the separation in the valence band closer to 150 me V— . 

The predicted strong increase in electron-hole Coulomb interaction in PbSe NWs should 
have major implications for other properties. This enhancement should increase the rate 
of the nonradiative Auger recombination as well as the rate of the inverse process, impact 
ionization. A high rate of impact ionization or efficient multiple exciton generation, combined 
with good conductivity that might be expected in PbSe NWs, suggests that these structures 
may be promising for photovoltaic applications. 

To summarize, we have developed a theory that describes both the energy spectra of 
individual electrons and holes and the absorption spectra of lead-salt NWs and NRs. Cal- 
culations show that even though spatial and dielectric confinement dramatically increase 
the exciton binding energy, the absorption spectra of PbSe NWs and NRs are practically 
unaffected, which should lead to insensitivity of these spectra to the surrounding media. The 
size dependence of lowest absorption peaks measured in PbSe NRs is very well described by 
the developed theory. It should be straightforward to apply this model to PbS and PbTe 
NRs. 
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Appendix A: Effect of anisotropy on the nanowire energy spectra 



The cylindrically symmetric Hamiltonian in Eq. fl4]) can be derived from the full Hamil- 
tonian in Eq. by transformation to the new coordinate system connected with NW 
direction. The full Hamiltonian is defined with respect to a crystallographic direction of the 
Brillouin zone, where the z-axis is pointed towards one of the L-points, and we will call 
this coordinate system the primed system, {x',y', z'} . We need to express Eq. ([3]) in the 
new coordinate system where the z-axis is directed along the rod axis, called the unprimed 
system, {x, y, z} . To do this, we use a coordinate rotation, and define the x-axis such that 
the rotation occurs in the x-z plane. In the rotation, vector quantities, such as p or & are 
transformed using the rotation matrix, p' = R(9)p , with R defined as 



R{9) 



cos 9 
1 



sin 9 




ysin^ cos 9 j 
This transformation expresses the squared momenta in Eq. 



(Al) 



as: 



~/2 
Px 

~/2 
Pz 



cos 2 ^-sin2^ 2 + sin 2 ^ 
sin 2 9p 2 x + sin 29p x p z + cos 2 9p 2 z 



(A2) 
(A3) 



and the diagonal and off-diagonal elements of the matrix of Hamiltonian in Eq. 
coordinate system as: 

1 



m t 



P'y 2 ) + -P^ 



P t a'J x + P t a' y p' y + Pia' z p' 



r' z p' z 



cos 2 9 sin 9 
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m t mi 

sin 2 9 cos 2 9 



P, 



1 -2 

—Py + 



mi 



P\ 



+ sin 20 ( — — ] p x p z 

\nii m t J 

(P t cos 2 9 + Pi sin 2 9)a x p x + P t cr y p y 
+(P t sin 2 9 + Pi cos 2 9)a z p z + 
+1- sin 29 (Pi - P t )(<y z Px + v x p z ) 



m new 



(A4) 



(A5) 



Notice that neither elements are cylindrically symmetric in the new coordinates. To en- 
force this symmetry, we rewrite these expressions in a form that separates a cylindrically 
symmetrical part, formally: aO x + bO y = (1/2) (a + b)(6 x + 6 y ) + (1/2) (a - b)(O x - 6 y ) . 
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The first term, which has cylindrical symmetry, is used in the zero-th order Hamiltonian, 
and the second term creates the asymmetric perturbation. This procedure produces the 
Hamiltonian in Eq. (J3J), along with the perturbation matrix 



/ 



x (| sin 2 6(pl - p 2 y ) + sin 26p x p z ) 

^(P~ Pt){ sin 2 9{& x p x - &yp y )-\- 
+ sin 29(a z p x + a x p z )} 



2 ^(P - P t ){sin 2 9(a x p x - a y p y ) + \ 
+ sin 26(a z p x + a x p z )} 



U U-4 



(A6) 



y + sin 26{a z p x + a x p z ) ) x (| sin 2 Q{p\ - p 2 y ) + sin 26p x p z ) j 

We study the effect of anisotropy described by Eq. (1A6j) on the energy spectrum of electrons 
and holes. Figure [10] compares the energy of the lowest electron levels in a 4 nm PbSe NW 
calculated within the cylindrical approximation and with complete numerical inclusion of 
the anisotropy. The anisotropy was taken into account by diagonalizing the matrix elements 
of H an in the space of the highest 20 valence and lowest 20 conduction states (that is, 
including the highest ten and lowest ten doubly degenerate electron and hole levels.) One 
can see in Fig. [10] that the anisotropy in PbSe splits the nearly degenerate energy levels, 
whose radial or angular quantum momentum numbers differ by one in radial or angular 
quantum momentum numbers, while necessarily leaving the Kramer's degeneracy unbroken. 
The splitting should broaden the energy levels without an overall shift in the level position. 



Appendix B: Calculations of the one dimensional Coulomb potential 

Calculation of the one dimensional Coulomb potential in Eq. ([TU]) and self interaction 
energy in Eq. f[T8|) can be greatly simplified by initial averaging over angular variables. For 
the U\ term of Eq. (|T9|) the angular integration results 

(U 1 )(z)= dp ePe dp h p h \® e \ 2 \$! h \ 2 Vi{p e , Ph ,z) , (Bl) 
Jo Jo 



e 2 ^ (z 2 + Pl + pl 



where 

Vi{p„p hi z) = -4vr— '—Q-i/2 ( - 9 re ' rh ) (B2) 

and Q n is the Legendre function of the second kind. The two remaining radial integrals are 
evaluated numerically. 
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FIG. 10. Effect of the energy spectrum anisotrpy on the energy of the ID subband bottom in a 4 nm 
PbSe NW grown along the (a) (111) and (b) (100) crystal axes. The "approximate" calculations 
are conducted within the cylindrical approximation, which gives Eq. (|15p for the energy levels. 
The "full" calculations are performed as described in the text. The energy levels are labeled by 
the angle between the L-point and the rod growth axis. Note that the 6 = energy levels do not 
require perturbation, as H an = for that angle. 

For the second term in Eq. (Tl9l) . XJi , the angular integrals vanish unless m = leaving 
only this term from the sum. This results in the following expression for (U2)(z) : 



ITT \( )- _8 — f°° A — - K ™) K o( Ru ) K l( Ru ) cos ( uz ) 

{ 2)[Z) ~ n K s J U K s h(Ru)K (Ru) + K m I (Ru)Ki(Ru) 
j dp e p e \ty e \ 2 I (up e )^ (J dp h p h \ty h \ 2 I (upi 



x 



(B3) 



i e (u) i h (u) 

To calculate the integrals i e and ih in Eq. flB3j) . we approximate the squared wavefunctions 
as a short sum of the form \^ e \ 2 = Yl n =i — Pe") 5 with N ^ 8 . Even with so few terms, 
the maximum relative error is typically < 10~ 7 . This allows us to solve these two integrals 
analytically: 

N 

iJu) 



y^An dp e p e (l - f^)I (up e ) 

n=l J( > 

Rh(u) R 2+2n xF 2 (1 + n; 1, 2 + n; R 2 u 2 /4) 



n=l 



U 



2 + 2n 



(B4) 



where p F q is the generalized hypergeometric function. The remaining integral over u in 
Eq. (IB3I) is performed numerically. 
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Lastly, the two self interaction terms in Eq. (ITS]) . U e and Uh , after angular integrations 
are reduced to 

2 e 2 f°° ( f R \ 
(U e ,h) = — z2 dui dp eyh p e , h \y e)h \ 2 I^(up e)h ) I X 

Ks m=0 Jo ^° / 

(k s - K m )K m (Ru)(K m _i(Ru) + K m+1 (Ru))(2 - £ m0 ) 

K s K m (Ru)(I m -i(Ru) + (Ru)) + {Ru){K m ^{Ru) + K m+1 (Ru)) ' v 7 

The two dimensional integrals in Eq. (1B5P was taken numerically. It is summed over only 
the first ~ 20 values of m , as the sum converges rapidly. 

Appendix C: Numerical calculation of the exciton binding in PbSe nanorods 

Our analytic model makes the assumption that the ID exciton is only weakly confined 
along the NR axis. In this case the finite length of the NR affects only the exciton center of 
mass motion. To verify this assumption, the ID Hamiltonian was numerically diagonalized, 
while treating both binding and confinement exactly. As an orthogonal basis for this diago- 
nalization we used a sufficiently large set of electron and hole plane waves that satisfied the 
single particle boundary conditions. The ID exciton wavefunction in this basis set can be 
written as: 

,Tr a 2 • fn e TTZ e \ . fn h TTZ h \ 

^ lD= 1^1^ A "«Af m {— — ) S111 {— — ) ( C1 ) 

n e =l n h =l 

where A ne ^ nh are the numerical coefficients. 

The kinetic energy is diagonal in this basis, and matrix elements of Eq. (|2T|) can be eval- 
uated analytically. Calculation time was dominated by evaluation of these matrix elements 
and scaled as 0( N e Nh ). For N e = Nh ~ 30 , calculations were sufficiently converged for the 
lowest few dozen states, and required roughly one minute of computation time on a desktop 
computer. 

Fig. [TTl shows the square of ID wavef unctions, |\l/iz)| 2 , calculated both the numerically 
and analytically as a function of z e and Zh- For the lowest two exciton states |\I/iz)| 2 shows 
good agreement between the numerical model and the analytical calculation. This is because 
the electron and hole are strongly localized around each other and do not feel the effects of 
confinement at the edges of box. As a result, the wavefunction orients along the coordinates 
associated with Coulomb binding, z and z cora , roughly along the graph diagonals. On the 
other hand, by the 17th excited state, also shown in Fig. [TTJ, the numerical and analytical 
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calculations disagree greatly. This is because the higher kinetic energy of this state causes 
the wavefunction to reach the edges of the box and feel confined. And, as a result, it begins 
to orient along the box coordinates, z e and Zh , associated with confinement. In general, our 
analytic model shows good agreement for the lowest ~ 10 states for each pair of nanowire 
bands. 



State #1 State #2 State #17 




electron position z c electron position z e electron position z c 



FIG. 11. Comparison of the numerically and analytically calculated ID exciton wavefunctions 
I^id) 2 . Each subplot has axes z e and ranging along the length of the nanorod from to L . 
The two lowest energy states and the 17th state are shown. 

Appendix D: Choice of the room temperature band parameters 

The absence of reliable room temperature energy band parameters for bulk PbSe has lead 
to several problems in the quantitative description of spherical PbSe NC electronic prop- 
erties within effective mass theory, and as a result, to some controversy on their electronic 
structure^^— . As has been noted^ 1 ^, effective mass theory significantly overestimates the 
energy gap in PbSe NCs (though not in PbS.) In addition, the nature of the 2nd optical tran- 
sition is still a source of debate^— , whether it is of symmetry type S-P or P-P. Considering 
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the body of experimental evidence, the explanation put forward by Franceschetti^ seems to 
offer the simplest explanation of this controversy, that the electron and hole P states are split 
into P± and P\\ states by the anisotropy of the bands, and the second transition is of type 
Pii - Pii . These two problematic aspects of experimental spectra of PbSe NCs for effective 
mass theory- overestimation of the bandgap and the symmetry of the 2nd transition- as 
well as the observation of several optical transitions in a wide range of energies can be used 
for extraction of a real set of the energy band parameters. 

Although the extraction of the set of energy band parameters from room temperature 
absorption spectra is possible, it is likely that many sets of parameters will equally well fit 
the first few optical transitions. In order to increase the accuracy of the fit, we want to 
somehow incorporate the energy band parameters in low temperature experiments in bulk 
PbSe. So, the total band edge effective masses for electrons and holes at T = 4 K are held 
constant at the values from experiment^. In addition, to limit the degrees of freedom in 
the fit, the anisotropy of the far-band contributions to both the electron and hole are held 
equal. That is, mf /mf = mj~ jm^ , even though their individual values will differ. With 
these constraints, a fit is performed using the body of literature data4^^~— for the first 
transition, and the data from Koole^ for the second and third transitions. 

The final set of room temperature parameters are shown in Table 
the set of low temperature parameters reported for bulk PbSe in Ref. 



V 



B] together with 



461 ] . The transition 



energies calculated using these parameters are shown in Fig. [12] The anisotropic effective 
mass calculations were performed using the method outlined in Ref. 64] and the results 
compared to the energies measured in Ref. 45[ , ignoring those points criticized in Ref. 44 ] 
as possibly being 2nd derivative artifacts. 
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FIG. 12. Calculations of the lowest electron levels in spherical PbSe NCs. (a) Splitting of the 
P state induced by the fully anisotropic Hamiltonian in a 4 nm radius NC. Anisotropic states 
are labeled by writing the state in the basis of isotropic states, and labeling it by the isotropic 
state with largest coefficient, (b) The size dependence of the transition energies in spherical PbSe 
NCs. Experimental data^i are shown by symbols. The solid lines show the size dependence of 
optically allowed transitions calculated in a fully anisotropic effective mass model. The optically 
allowed transitions occur between the states of the same symmetry but opposite parity, and we 
label them by a symmetry type, which is common for both states. Open points indicate transitions 
originating from the L-point in the Brillouin zone, while half-open points are suggested to be 



from the S point as in Ref. [45(. The dashed line shows the size dependence of lowest confined 
level connected with the £ point of the Brillouin zone, calculated in a parabolic effective mass 
approximation as explained in the text. 
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